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Abstract 

In this paper we present an algorithm for enumerating without repetitions all the non- 
crossing generically minimally rigid bar-and-joint frameworks under edge constraints (also 
called constrained non-crossing Laman frameworks) on a given generic set of n points. 
Our algorithm is based on the reverse search paradigm of Avis and Fukuda. It generates 
each output graph in 0(n ) time and 0(n) space, or, slightly different implementation, in 
0(n 3 ) time and 0(n 2 ) space. In particular, we obtain that the set of all the constrained 
non-crossing Laman frameworks on a given point set is connected by flips which restore 
the Laman property. 

Key words: geometric enumeration; rigidity; constrained non-crossing minimally rigid 
frameworks; constrained Delaunay triangulation. 



1 Introduction 

Let G be a graph with vertices {1, . . . ,n} and m edges. G is a minimally rigid graph (also 
called Laman graph) if m = 2n — 3 and every subset of n' < n vertices spans at most 2n' — 3 
edges. An embedding of the graph G(P) on a set of points P = {pi, • • • ,p n } C R 2 is a 
mapping of the vertices to points in the Euclidean plane i i— > pt 6 P. The edges ij of G are 
mapped to straight line segments PiPj- An embedding is planar or non-crossing if no pair 
of segments PiPj and pkPi corresponding to non-adjacent edges i,j {k,l} have a point in 
common. 

An embedded Laman graph on the generic point sets is called Laman framework and has 
the special property of being minimally rigid [16,21], when viewed as bar-and-joint frame- 
works with fixed edge-lengths, which motivates the tremendous interest in their properties. 
Let F be a set of non-crossing edge (bar) on P. A Laman framework containing F is called 
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F- constrained. In this paper we give an algorithm for enumerating all the F- constrained 
non-crossing Laman frameworks. 

In this paper we only consider planar embedded Laman graphs, namely Laman frame- 
works. Then we simply denote a vertex pi by i and an edge piPj by ij. 

Novelty. To the best of our knowledge, this is the first algorithm proposed for enumer- 
ating (without repetitions, in polynomial time and without using additional space) all the 
-F-constrained non-crossing Laman frameworks. We achieve 0(n 3 ) time per graph in 0(n 2 ) 
space by using reverse search. For the unconstrained case, using a different method, the 
same time and space complexity results we obtained in [5]. The method presented here uses 
matroid property and is considerably simpler than [5]. 

Historical Perspective. The reverse search enumeration technique of Avis and Fukuda 
[3, 4] has been successfully applied to a variety of combinatorial and geometric enumeration 
problems. The necessary ingredients to use the method are an implicitly described connected 
graph on the objects to be generated, and an implicitly defined spanning tree in this graph. 
In this paper we supply these ingredients for the problem of generating Laman frameworks. 

Relevant to the historical context of our work are the results of Bereg [8, 11] using reverse 
search combined with data-specific lexicographic orderings to enumerate triangulations and 
pointed pseudo-triangulations of a given point set. We notice in passing that there exist sev- 
eral other algorithms for enumerating (pseudo-)triangulations [1, 12, 14], but they are based 
on different techniques. 

Also relevant is the pebble game algorithm of Jacobs and Hendrickson [17] for 2-dimensional 
rigidity, see also [9]. Our complexity analysis relies on the recent results, due to Lee, Streinu 
and Theran [22,23], regarding the detailed data-structure complexity of finding and main- 
taining rigid components during the pebble game algorithm. Indeed, the time-space trade-off 
of our algorithm is inherited from [23]. 

Related previous works. To better put our problem (and solution) in context and 
relate it to previous work, we briefly discuss now the difference between generating arbitrary 
Laman frameworks, as opposed to just pointed pseudo-triangulations. 

A pointed pseudo-triangulation is a special case of a non-crossing Laman framework 
[27], where every vertex in the embedding is incident to an angle larger than ir. Pseudo- 
triangulations are connected via simple flips, in which the removal of any non-convex-hull 
edge leads to the choice of a unique other edge that can replace it, in order to restore the 
pointed pseudo-triangulation property. The flip graph of all pointed pseudo-triangulations is 
a connected subgraph of the graph of all the Laman frameworks. In fact, it is the one-skeleton 
of a polytope [25], and the reverse search technique can be directly applied to it. Bereg's 
efficient algorithm makes use of specific properties of pointed pseudo-triangulations which do 
not extend to arbitrary non-crossing Laman frameworks. In particular, remove-add flips are 
not unique, relative to the removed edge, in the case of a non-crossing Laman framework. 

In our previous paper [5] we showed that a set of all non-crossing Laman frameworks on 
a point set is connected via these flips. The proof relies on some properties of one-degree-of- 
freedom mechanisms which is obtained when removing one edge from the Laman frameworks, 
which cannot extend to edge constrained case. 

In this paper we use triangulation as an important tool. We use the fact that (constrained) 
non-crossing Laman frameworks form the bases of a matroid defined on any triangulation of P 
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to define a parent function. The edge constrained case appears in an application to structural 
engineering as will be seen below, and the proposed algorithm will much reduce the output 
size compared with the unconstrained case. Therefore, we can enumerate non-crossing Laman 
frameworks for previously unsolved real problems. 

Motivating application. We describe now briefly how this problem came to our at- 
tention via the work of the third author. Graph theoretical approaches are widely used in 
structural mechanics [19], where the edges and vertices in the graph represent the bars and 
rotation- free joints of a structure called a truss. It is well-known that the stiffest truss under 
static loads is statically determinate that is equivalent to a Laman graph [7]. 

Another bar-joint system, which is widely used in industrial application, is a link mecha- 
nism that is unstable and generates large deformation or changes the direction of the nodal 
displacement. The applications of link mechanisms can be found in, e.g., suspension of auto- 
mobiles, robot hands, umbrellas, crank shafts, etc. Kawamoto et al.[20] presented a method 
by enumeration of the planar graphs to find an optimal mechanism. However, their method 
is developed for their specific problem, and no general approach was given. 

Recently, a new type of mechanism called compliant mechanism has been developed and 
applied mainly in the field of micro-mechanics. A compliant mechanism has flexible parts to 
stabilize the structure, which is contrary to the conventional unstable mechanism. Although 
a compliant mechanism is usually modelled as a continuum with elastic joints, it is possible to 
generate the similar mechanism by a bar-joint system. Ohsaki and Nishiwaki [24] presented a 
method for generating flexible bar-joint mechanisms using nonlinear programming approach, 
and found that the optimal structure is statically determinate, i.e. minimally rigid. They 
utilized snapthrough behaviour to generate multi-stable mechanism that has multiple self- 
equilibrium states. Such mechanism can be used for switching device, robot hand, gripper, 
deployable structure, etc. In their method, the optimal locations of bars and joints are found 
from the highly connected initial structure that has bars between all the pairs of the nodes 
whose distances are small enough. Combining an implementation of our our earlier method 
for generating unconstrained Laman frameworks [5] with this nonlinear programming ap- 
proach, we could obtain many new compliant mechanisms with up to 10 joints [18]. However 
the number of Laman frameworks grows too rapidly to allow a complete enumeration for 
significantly larger examples. 

In view of practical requirements, the optimal structure should not have intersecting mem- 
bers, and some pre-selected members should always exist. Therefore, the computational cost 
may be much reduced if the candidate set of statically determinate non-crossing trusses (non- 
crossing Laman frameworks) are first enumerated. Thus, it is desirable to enumerate Laman 
frameworks for which the existing edges are specified, and this is the problem considered in 
the current paper. 

2 Preliminaries 

Let L be a non-crossing Laman framework on a given point set P. A mechanism is a flexible 
framework obtained by removing one or more edges from a Laman framework. Its number of 
degrees of freedom or dofs, is the number of removed edges. We will encounter mostly one- 
degree- of -freedom (ldof) mechanisms, which arise from a Laman framework by the removal 
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of one edge. In particular, a mechanism with k dofs has exactly 2n — 3 — k edges, and each 
subset of n' vertices spans at most 2n' — 3 edges. A subset of some n' vertices spanning 
exactly 2n' — 3 edges is called a rigid block. A maximal block is called a rigid component or a 
body. 

The Laman frameworks on the generic point set with n points form the set of bases of 
the generic rigidity matroid (or simply called Laman matroid) on K n , see [16]. The bases 
have all the same size 2n — 3. Bases may be related via the base exchange operation, which 
we will call a flip between two Laman frameworks. Two Laman frameworks L\ and L2 are 
connected by a flip if their edge sets agree on 2n — 4 elements. The flip is given by the pair of 
edges (ei, e-i) not common to the two bases, e± £ L\ — L2, e<i € L2 — L\. Using flips, we can 
define a graph whose nodes are all the Laman frameworks on n vertices, and whose edges 
correspond to flips. 

It is well-known that the graph whose nodes are the bases of a matroid connected via 
flips, is connected. But a priori, the subset of (constrained) non-crossing Laman frameworks 
may not necessarily be. We will prove this later in Section 3. 

Reverse search is a memory efficient method for visiting all the nodes of a connected graph 
that can be defined implicitly by an adjacency oracle. It can be used whenever a spanning 
tree of the graph can be defined implicitly by a parent function. This function is defined for 
each vertex of the graph except a pre-specified root. Iterating the parent function leads to a 
path to the root from any other vertex in the graph. The set of such paths defines a spanning 
tree, known as the search tree. 

3 Constrained Non-crossing Laman Frameworks 

Let T be a triangulation on a given set of n points P in the plane containing k triangles. The 
angle vector of T is the vector of 3k interior angles sorted into non-decreasing order. Let F 
be a non-crossing edge set on P. An F-constrained triangulation on P is one that contains 
F as a subset. Many facts about these triangulations are contained in the survey by Bern 
and Eppstein [10]. If F is an independent set in the Laman matroid on K n , then a Laman 
framework containing F is called F-constrained. The following lemma clearly holds from the 
known fact about matroids (see, e.g. [28]). 

Lemma 1. Let F be a non-crossing edge set on P that is an independent set in the Laman 
matroid on K n . Let us fix F-constrained triangulation T on P. Then a set of bases of the 
Laman matroid on T whose edge set contains F forms matroid. 

Lemma 2. Let F be a non-crossing edge set on P that is an independent set in the Laman 
matroid on K n . Every F-constrained triangulation T on P contains an F-constrained Laman 
framework. 

Proof. Since T is statically rigid, it is generically rigid, and hence contains a base B of the 
Laman matroid (see, e.g. Whiteley [29]). Since F is independent in the Laman matroid, it 
can be extended to a base by adding edges from B — F. □ 

Two points a and b are visible if no edge of F properly intersects the segment ah. ah is 
visible to point c if the triangle abc is not properly intersected by an edge of F. 
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Definition 1. (CDT) An F-constrained Delaunay Triangulation (CDT) contains the edge 
ab between points a and b in P if and only if a is visible to b, and some circle through a and 
b contains no point of P visible to segment ab. We call ab a D(elaunay)-edge. (Definition 1, 
[10]) 

If P has four or more co-circular points, using a linear transformation as described in 
[6], we may transform P into a point set P with a unique CDT. P and P have the same 
non-crossing Laman frameworks. Then we will assume in what follows that P has a unique 
CDT. 

Definition 2. (D-flip) Let ac be an edge ofT — F which is the diagonal of a convex quadri- 
lateral abed contained in T . The edge flip which replaces ac by edge bd is a D-flip if the 
circumcircle of the triangle abc contains the point d in its interior. ( Equivalently the circum- 
circle of the triangle acd contains the point b in its interior.) We call ac a F -illegal edge, 
otherwise it is called F -legal edge. 

Note that a D-flip increases the angle vector lexicographically. This can be used to prove 
the following. 

Fact 1. The CDT has the lexicographically maximum angle vector of all F-constrained tri- 
angulations on P. (Theorem 1, [10]). 

Fact 2. An F-constrained triangulation T can be converted to the CDT by at most 0{n 2 ) 
D- flips, taken in any order. (Lemma 4, [10]). 

Now let us consider non-crossing Laman frameworks. 

Definition 3. A Constrained Delaunay Laman Framework ( CDLF) is an F-constrained 
Laman framework that is a subset of the CDT. 

Definition 4. (L-flip) An L-flip is an edge insertion and deletion that takes a Laman frame- 
work L to a new Laman framework V . 

Theorem 1. Every F-constrained non-crossing Laman framework L can be transformed to 
a CDLF by at most 0(n 2 ) L-flips. 

Before giving the proof of Theorem 1, we define the underlying triangulation T(L) for a 
non-crossing Laman framework L, which is L-constrained Delaunay triangulation. 

Definition 5. Let L be a non-crossing Laman framework. The L-constrained Delaunay 
triangulation T{L) is constructed by adding edges to L in the following way. First, we add 
convex hull edges if missing in L, and for each planar face of L we compute its (internal) 
constrained Delaunay triangulation, and add these new edges. 

By this construction, we have the following fact: 

Fact 3. Let L be a F-constrained non-crossing Laman framework. Then all F -illegal edges 
in T(L) are the edges of L — F . 
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Proof of Theorem 1. Suppose that L has an F-illegal edge. Then, from Fact 3, such a F- 
illegal edge, say ac, is among the edges of L — F. Consider removing ac. The updated 
underlying triangulation T(L — ac) contains L — ac, and by Lemma 2 and the matroid property 
shown by Lemma 1, there is some edge st in T(L — ac) such that V = L — ac + st is a non- 
crossing Laman framework. The fact that ac is F-illegal edge implies that there exists at 
least one D-fiip when updating T(L) to T(L — ac), which lexicographically increases the angle 
vector. From Fact 1 and 2, repeating this procedure 0(n 2 ) times, we eventually reach the 
F-constrained Delaunay triangulation T(L"), and L" is the required CDLF. □ 

For edges e = ij with i < j and e' = kl with k < I, we use the notation e -< e' or 
e' >- e when e is lexicographically smaller than e' i.e., either i < k or i = k and j < I, 
and e = e' when they coincide. For an edge set A we use the notations max{e | e G i} 
and min{e | e E A} to denote the lexicographically largest and smallest labelled edges in A, 
respectively. 

Definition 6. (Lexicographic edge list) Let E = {e± -< ei -< . . . -< e m } and E' = {e'i ~< e' 2 -< 
. . . -< e' m } be the lexicographically ordered edge lists. Then E is lexicographically smaller than 
E' if ej -< e[ for the smallest i such that e.{ ^ e[ . 

Theorem 2. The set of F- constrained non-crossing Laman Frameworks on a point set P is 
connected by 0(n 2 ) edge flip operations. 

Proof. By the previous theorem, from any non-crossing Laman framework we can perform 
L-flips 0(n 2 ) time to reach a CDLF, say L. Let L* be the CDLF with lexicographically 
smallest edge list. We show that we can do edge flips from L to L* , at most n — 3 times, 
maintaining the non-crossing Laman property. Indeed, delete from L the largest indexed 
edge ac in L — L* . By the matroid properties, there will always be an edge in L* — L to 
insert that maintains the Laman framework property. In fact, the ldof mechanism L — ac 
is not maximum component and can be extended to a base by adding an edge from L* — L. 
Planarity is maintained since both L and L* are subsets of the CDT which is non-crossing. 

Let n e be the number of edges of the triangulation on P. By Euler's formula, n e = 
3n — h — 3<3n — 6 holds, where h > 3 is the number of points on the convex hull of P. 
Then L has at most n — 3 edges which are not in L*. Since we can replace one of such edges 
by an edge of L* by one edge flip, we eventually reach L* by n — 3 edge flips. □ 

4 Algorithm 

Let £ be a set of F-constrained non-crossing Laman frameworks on P, and T>£ be a set of 
the CDLFs. Clearly VC C £ holds. Let L* be the CDLF with lexicographically smallest 
edge list. We define the following parent function f : C — {L*} — » C based on Theorem 1 and 
2. 

Definition 7. (Parent function) Let L G C with L ^ L* . L' = L — ac + st is the parent of 

L, where 

Case 1: L £ VC, 

ac = max{e | e € L — L*} and st = min{e €L*— L\L — ac + e £L £}, 
Case 2: L£ C- VC, 



6 





Figure 1: The example of the parent function. The solid lines represent edges of L, and the 
dotted lines represent additional edges for the underlying triangulations. 



ac = max{e G T(L) \ e is F -illegal edge in T(L)} and st = min{e G T(L — ac) — (L — ac) 
L — ac + e G £}. 

To simplify the notations, we denote the parent function depending on Case 1 and Case 
2 by fi : T>L — {L*} — > DC and / 2 : C — T>C — > C, respectively. In Fig.l, we show the 
example of / 2 (£) in which L is not a CDLF: removing the largest indexed F-illegal edge 
37, and inserting the smallest indexed edge 12 in T(L — ac) — (L — ac), we obtain another 
non-crossing Laman framework shown in the rightmost and upper corner of Fig. 1. 

The reverse search algorithm can be considered on the underlying graph in which each 
vertex corresponds to a non-crossing Laman framework and two frameworks are adjacent if 
and only if one can be obtained from the other by a L-fhp. Then, for L' G C the local search 
is given by an adjacency function, Adj, defined as follows: 



Adj(L',ei,e 2 ) - 



V — ei + e 2 if L' - e\ + e 2 G £, 
null otherwise, 



where e\ G L' — F and e 2 G K n — L' . The number of candidate edge pairs (ei, e 2 ) is 0(n 3 ). 

Let elist^' and elist/^ n be the list of edges of L' and K n ordered lexicographically, let S(L') 
and 5{K n ) be the number of elements of elist^' and elist/^ n and let elist^/(i) and elist^ n (i) 
be the i-th elements of elistz/ and elist^: n , respectively. We also denote the above defined 
adjacency function by Adj(L' ,i,j) for which e\ = elist^(i) with e\ F and e 2 = e\'\stK n (j) 
with e 2 ^ L'. Then, based on the algorithm in [3, 4], we describe our algorithm in Fig. 2. 

As we will show later, both the parent function and the adjacency function need 0(n 2 ) 
time for each process. Then, the while-loop from line 4 to 17 has S(L') ■ 5(K n ) iterations 
which require 0(n 5 ) time if simply checking the line 8 and 10. In order to improve 0(n 5 ) 
time to 0(n 3 ) time we claim the following two lemmas: 

Lemma 3. Let L and L' be two distinct CDLFs which are subgraphs of T(L*) for which 
L = Adj(L' , ei, e 2 ) for e\ G V — F and e 2 G K n — L' . Then, fi(L) = L' holds if and only if 
e\ and e 2 satisfy the following conditions: 
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Algorithm Enumerating F-constrained non-crossing Laman frameworks. 

1: L* := the CDLF with lexicographically smallest edge list; 
2: L' := L*; i,j := 0; Output(L'); 
3: repeat 

4: while i < S(L') do 
5: i := i + 1; 
6: while j < 5{K n ) do 
7: j := j + 1; 

8: if e\\st L /(i) £ F, elistft- n (j) ^ L' and Adj(L',i,j) / null then 

9: L := Adj(L',i,j); 

10: if /i(L) = L' or f 2 (L) = L' then 

11: L' := L; i,j := 0; 

12: Output(L'); 
13: go to line 4; 

14: end if 

15: end if 

16: end while 
17: end while 
18: if L' ^ L* then 
19: L := L'; 
20: if L € VC then L' := 
21: else L' := / 2 (^); 

22: determine integers pair (i, j) such that Adj(L' = L; 
23: i:=i — l; 
24: end if 

25: until V = L*,i = S(L') and j = 8(K n ); 

Figure 2: Algorithm for enumerating F-constrained non-crossing Laman frameworks. 

(a) ei G L*, 

^ e 2 G T(L*) - L* ; 

(c) ei ^ min{e G L* - L' \ L' - a + e G C}, 

(d) e 2 J- max{e \ e £ L' - L*}. 

Lemma 4. Lei L and L' be two distinct non-crossing Laman frameworks for which L = 
Adj(L' , e±, e2) for edges e\ G L' — F and e 2 G K n — V , and let L G C — DC. Then, f2(L) = L' 
holds if and only if e\ and e 2 satisfy the following conditions: 

(a) e\ is F -legal edge in T(L'), 

(b) e 2 €K n -T(L'), 

(c) ei -< min{e G T(L') - V \ L' - e x + e G £}. 

(d) e 2 = max{e G T(JJ — e\ + e 2 ) \ e is F -illegal edge in T(L' — e\ + e 2 )}. 
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We will explain later (in the proof of Theorem 3) how Lemmas 3 and 4 are used to obtain 
0(n 3 ) time for generating each output of our algorithm. Notice that for V and L G C such 
that L = L' — ei + e2, at most one of f\(L) = L' and / 2 (L) = V holds from the conditions 
(b) of Lemmas 3 and 4. 

Proof of Lemma 3. ("only if'-part.) Since fi(L) = L' holds, e\ and e 2 must be chosen as st 
and ac in Case 1 of Definition 7. From Definition 7, e 2 (= ac) £ L — L* holds. Since L G VC, 
L C T(L*) and e 2 G T(L*) - L* holds, thus we have (b). Similarly since ei(= st) £ L* — L C 
L* , we have (a). 

From ei = st, we have 

L' — e\ = (L — ac + st) — e\ = L — ac. (1) 

Let e' = min{e £ L* — L' \ V — e\ + e £ C}. Suppose (c) does not hold. Then e' -< e\ 
holds. (Note that the equality does not hold since c\ € V — F.) Then from Eq.(l) and 
e' -< e\ = st -< ac (which comes from Definition 7), 

e' = min{e G L* - L' \ V - a + e G £} 

= min{e £ L* — (L — ac + st) \ L — ac + e G £} (from (1)) 

= min{e G L* — L \ L — ac + e G £} (from e' -< st -< ac). 

Thus, e' would have been selected instead of c\ when the parent function f\ is applied to L, 
which contradicts e\ = st. Hence, (c) holds. 

Finally, let e" = max{e | e G L' — L*}, and suppose that (d) does not hold and e 2 -< e" 
holds. (Note that the equality does not hold since e 2 ^ L'.) Since st -< ac = e 2 -< e", 

e" = max{e \ e €. L' — L*} 

= max{e | e G (L — ac + st) — L*} 
= max{e | e G L — L*}. 

Then e" would have been selected instead of e 2 when the parent function f\ is applied to L, 
which contradicts e 2 = ac. Thus, (d) holds. 

("if'-part.) From (a) and (b), L = L' — e± + e 2 G VC. Since e\ G L* from (a), (d) implies 

e 2 >- max{e | e G L' — L*} 

= max{e | e G (L + e\ — e 2 ) — L*} 
= max{e | e G (L — e 2 ) — L*}. 

Thus, e 2 = max{e \ e £ L — L*} holds, and hence f\ chooses e 2 for an edge ac to be deleted 
from L. 

From this we have L — ac = L' — e\ + e 2 — ac = L' — e\ . Since e 2 ^ L* from (b) , (c) implies 

ei -<; min{e G L* - L' | L' - e\ + e G £} 

= min{e £ L* — (L + e\ — e 2 ) | L — ac + e G £} 
= min{e £ L* — (L + e\) \ L — ac + e £ £}. 

Since e\ £ L* — L, e\ = min{e £ L* — L \ L — ac + e £ C}. Thus, f\ chooses e\ for an edge 
to be added, and fi(L) returns V . □ 
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Proof of Lemma 4- ("only if'-part.) Since f2(L) = L' holds, e\ and e2 must be chosen as st 
and ac of Definition 7. From e2 = ac, we have 

L — ac = L' — e\ + e<i — ac = L' — e±. (2) 

Since st G T(L — ac) — (L — ac) from Definition 7 and Fact 3, st is F-legal in T(L — ac) and 
hence F-legal in T(L — ac+ st)(= T(L')). Thus, from e\ = st, (a) holds. Since e\ is F-legal 
in T(L'), we have 

T(L') = T{L' -a) =T{L-ac). (3) 

Let e' = min{e G T(L') — L' \ L' — c\ + e G £}. Suppose that (c) does not hold. Then e' -< ei 
holds. (Note that the equality does not hold since e\ E L' — F.) From Eq.(2) and (3), 

e' = min{e G T(L') - L 1 | L' - e x + e G £} 

= min{e G T(L — ac) — (L — ac + e±) \ L — ac + e G £} 
= min{e G T(L — ac) — (L — ac) \ L — ac + e G £}. 

Then, e' would have been selected when the parent function is applied to L, which contradicts 
c\ = st. Hence (c) holds. 

Next let us consider e2- Since e2(= ac) is F-illegal in T(L) from Definition 7, we have 
e 2 i T(L - ac) = T(L'). Thus (b) holds. 

Finally, (d) must hold since the parent function remove an lexicographically largest illegal 
edge in T(L). 

("if'-part.) From (a), c\ is F-legal in T(L'). Then, we have T(L') = T(L' — ei). Suppose 
that e2 is F-legal in T(L) = T(L' — c\ + e-i), T(L' — e\ + e2) = T(JJ — c\) holds and e2 is 
in T(V — ei) = T(L'), which contradicts the condition (b). Therefore, C2 is F-illegal edge in 
T(L), and (d) says that C2 is the lexicographically largest F-illegal edge in T(L). Thus, fi 
chooses C2 for an edge ac to be deleted from L, and L' — c\ = L — ac holds. 

From V — e\ = L — ac and the condition (c), e\ -< min{e G T(L — ac) — (L — ac + e\) \ 
L — ac + e G £} holds. Thus, e\ = min{e G T(L — ac) — (L — ac) \ L — ac + e G £}. (Note 
that ei G T(L - ac), because T(L - ac) = T(L' - a) = T(L') and a G T(L') hold from (a).) 
Thus, /2 chooses c\ for an edge to be added, and f2{L) returns V . □ 

By Lemmas 3 and 4, we will describe an 0(n 3 ) algorithm in the proof of the following 
Theorem 3. Before it, we give you the simple observation for checking the condition (d) in 
the Lemma 4 efficiently: 

Observation 1. Let T(L') be a L' -constrained Delaunay triangulation constrained by edges 
of a non-crossing Laman framework L' , and let e\ G L' be a F-legal edge in T(L') and 
?2 G K n — L' be an edge intersecting no edge of L' . Then T(L' — e\ + e2) = T(L' + e2) holds. 

Proof. This comes from the following two facts that: 

(1) T(L') = T(L' - ei) holds, and 

(2) ei is F-legal in T(L' + e 2 ). 
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The fact (1) clearly hold since e± is F-legal edge in T(L'). Using Lemma 4.2 in [15] the 
fact (2) immediately follows. Let us show how to update T(L') to T{L' + e 2 ). The edge e 2 
intersects some edges in T(L') — L', and we denote a set of such edges by C e2 . First, we 
delete the edges in C e2 . The resulting graph T(L') — C e2 has a hole bounded by the polygon 
Q. When inserting e 2 into T(L') — C e2 , e 2 splits Q into two polygon Qi and Q2- Then, 
T(L' + e 2 ) = (T(L') - C e2 ) + T(Qi) + T(Q 2 ) holds from the fact in [15], i.e. T(L' + e 2 ) is 
obtained from T(L') by (Delaunay) triangulating only Q\ and O2 independently. This implies 
that e\ never flips even if e\ is not a constrained edge. (Note that e 2 does not intersect ei 
from the lemma assumption.) Hence, e\ is F-legal in T(L' + 62). □ 

Theorem 3. The set of all F -constrained non-crossing Laman frameworks on a given point 
set can be reported in 0(n 3 ) time per output using 0{n 2 ) space, or 0(n 4 ) time using 0(n) 
space. 

Proof. As described in Section 3, we use a linear transformation if necessary to get a unique 
CDT. The complexity of testing the uniqueness of a CDT is 0(n 2 ) by simply testing the 
circumcircle of each triangle in the CDT to see there is another point other than vertices of 
the triangle on the circumcircle. 

Given a non-crossing Laman framework L' G C and L'-constrained Delaunay triangulation 
T(L'). The algorithm will check f^AdftL', e 1 ,e 2 )) = V or f 2 (Adj(L f , e 1 ,e 2 )) = V at line 10 
depending on the edge pair (ei, e 2 ). Here we will show that each condition in Lemmas 3 and 
4 can be checked in O(l) time for each of the 0(n 3 ) edge pairs (e\,e 2 ) by the following way. 

First, for all edges e 2 G elist^, we calculate the number of edges e\ G L' intersecting 
e 2 , which we denote by cross_n(e2, L'). If cross_n(e2, L') > 1, we delete e 2 from elist^- n since 
V — ei + e 2 is never non-crossing for any e\ G elist^/. If cross_n(e2, L') = 1, we associate a 
pointer of the edge e\ intersecting e 2 with e 2 , and we denote such e\ by cross_e(e2, L'). 

Next, for each e\ G elist/,/, we attach two flags to e\ which represent that e\ satisfies the 
condition (a) of Lemmas 3 and 4, respectively. This preprocessing can be done in 0{n) time 
by simply scanning all edges of elist^y. By this we can check the condition (a) in Lemmas 
3 and 4 in 0(1) time. Similarly, we attach two flags to e 2 G elist^ n which represent that e 2 
satisfies the condition (b) of Lemmas 3 and 4. This process can be done in 0(n 2 ) time. By 
this we can check the condition (b) in Lemmas 3 and 4 in 0(1) time. We can calculate the 
lexicographically largest edge in V — L* in O(n) time. By this we can check the condition 
(d) in Lemmas 3 in 0(1) time. 

Now let us consider how to identify a set of edges e 2 G elist^/ satisfying the condition (d) in 
Lemma 4 in 0(n 3 ) time with 0(n 2 ) space. (In the case of 0(n 4 ) time algorithm this process 
must be skipped, and the condition (d) in Lemma 4 will be checked simply by updating 
T(L') to T(L' — e\ + e 2 ) for each pair (ei,e 2 ) using O(n) time and O(n) space by applying 
the algorithm by Chin and Wang[13].) It can be done regardless of the removing edge e\ if the 
condition (a) in Lemma 4 holds. From Observation 1, we can say that the condition (d) holds 
if and only if e 2 is lexicographically largest F-illegal edge in T{L + e 2 ) when cross_n(e2, L') = 
0. And it is sufficient to check the condition (d) only in T(L— cross_e(e2, L') + e 2 ) when 
cross_n(e2, L') = 1. Updating the Delaunay triangulation takes O(n) time (see [2, 13, 15] for a 
linear time update of the constrained Delaunay triangulation). Thus we can attach a flag to 
e 2 G elistft- n in O(n) time which represents whether e 2 satisfies the condition (d) of Lemma 4 
or not, thus this preprocessing for all edges in elistx n takes 0(n 3 ) time. 
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By using the above mentioned data, we will show that for a fixed e\ G elistjy, the inner 
while-loop from line 6 to 16 can be executed in 0(n 2 ). In order to efficiently test the condition 
(c) of Lemmas 3 and 4, we prepare the data structure proposed by Lee, Streinu and Theran 
[22,23] in 0(n 2 ) time for maintaining rigid components obtained by deleting e±. This data 
structure supports a pair-find query which determines whether two vertices are spanned by 
a common component in 0(1) time using 0(n 2 ) space, or 0(n) time with O(n) space. From 
this, we can calculate Adj(V , e±, e 2 ) (i.e., determine whether V — e\ + e 2 G C) in 0(1) time 
with 0{n 2 ) space, or O(n) time with O(n) space, for each edge e 2 G elist^ n . Also, we can 
compute e' = min{e <G L* -V \ L'-ei + e G £} and e" = min{e G T(L')-L' \ L'-ei + e G £} 
in 0(n) time with 0(n 2 ) space, or 0(n 2 ) time with O(n) space. From e' and e" we can check 
condition (c) in Lemmas 3 and 4 in 0(1) time. 

We have to determine which of f\(Adj(L', e\, e 2 )) = L' and f2{Adj(L', e±, e2)) = L' should 
be checked for every edge e 2 G elistx„- This is simply done by checking whether e 2 G T(L*) — 
L* or e 2 e K n - T(L'). When e 2 G T(L*) - L*, we check whether f^AdjiL', ei, e 2 )) = L' 
holds or not. When e 2 G K n — T(L'), we check whether f2(Adj(L' , e±, e 2 )) = L' holds or not. 
Both can be done in O(l) time with On 2 space, or O(n) time with 0(n) space. 

By using the above mentioned data structure for maintaining the rigid components, we 
can perform both parent function and adjacency function in 0(n 2 ) time with 0(n 2 ) space, 
or 0(n 3 ) time with O(n) space. Thus, we have an 0(n 3 ) algorithm using 0(n 2 ) space, or 
0(n 4 ) algorithm using 0(n) space. □ 

Figure 3 illustrates an example of the search tree of a set of F-constrained non-crossing 
Laman frameworks on six points, where F = {13, 15, 26, 45, 56}. 

5 Conclusions 

We have presented an algorithm for enumerating all the constrained non-crossing Laman 
frameworks. We note in passing that the techniques in this paper can also be used to generate 
all F-constrained non-crossing spanning trees of a point set since they also form bases of the 
graphic matroid on any triangulation of P. The unconstrained case was considered in [4]. 
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